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Abstract 

We show a method to solve the problem of the brachistochrone as 
well as other variational problems with the help of the soap films that 
are formed between two suitable surfaces. We also show the interesting 
connection between some variational problems of dynamics, statics, optics, 
and elasticity. 



1 Introduction 

The calculus of variations can be used in the formulation of most of the physical 
problems. Its origin was the famous problem of the brachistochrone, the curve of 
shortest descent time.^ Johann BernouUi found that the curve solution to this 
problem is a cycloid, that is, the curve described by a point on the circumference 
of a circle that rolls along a straight line.^ 

The calculus of variations was also present in the solution of other classical 
problems, such as the catenary and the isoperimetric problem. Its development 
was parallel to that of mechanics, geometrical optics, and elasticity. In this 
article we show how these problems can be related with the Plateau problem, 
that is, the problem of finding the form that a soap film adopts for certain 
contour restrictions. "^^^ We also show how to design simple experiments with 
soap films to find the solution to other variational problems such as a generalized 
Steiner problem or the problem of the quickest approach. 

This article was inspired by reading the book "Demonstrating Science with 
Soap Films" by D. Lovett, in which a similar argument is used to generate 
the catenary.^ The paper can be useful for students and teachers of Physics 
undergraduate courses. A complete study of this subject requires the use of 
the Euler-Lagrange equations. However, given the intuitive character of the 
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Figure 1: The velocity of a heavy body after descending an altitude y along 7 
is y/2gy. 




analogies between mechanics, optics, and elasticity problems that we establish, 
together with the easy visualization of the solutions with soap films, this article 
may also be useful for students without knowledge of differential equations. In 
the Appendix we propose a problem to use the concepts that appear in the 
article in a new situation. 

2 Relation between the Brachistochrone and 
Plateau's problems. 

Johann Bernoulli posed the problem of the brachistochrone as follows (see Fig- 
ure 1): 

Let two points A and B be given in a vertical plane. To find the curve that a 
point M, moving on a path AMB, must follow that, starting from A, it reaches 
B in the shortest time under its own gravity} 

Denote by 7 a curve joining A and _B, and let us approximate this curve by 
a polygonal A = Pi < ■ ■ ■ < < Pj+i < • • • < P„ = B, such that the distance 
|PiPi+i| is very small, thus we can assume that the velocity along PiPj+i is 
constant and can be taken as the velocity of a heavy body after descending an 
altitude yi, i.e. Vi = \/2gyi, where g is the gravity acceleration and yi is the 
ordinate of the midpoint between P, and P^+i (see Figure 1). 

Then the total time of descent along the polygonal is given by: 
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The total time of descent along 7, t{'^), is then obtained from this approximation, 
taking the limit as the length of the segments of the polygonal goes to 0, thus 



Figure 2: Vertical surface E that aproximates the soap film bounded by the 
surfaces z = 0, z = and two vertical segments. Its orthogonal projection is 
the cycloid 7. 



The problem is to find the curve 7 that minimizes the above expression. The 
solution to this problem is given by solving a differential equation, known as 
Euler- Lagrange equation, whose solution is a cycloid.^ 

We are going to study a Plateau problem that, from a mathematical point of 
view, is equivalent to the Brachistochrone problem and, therefore, its solution 
is also the cycloid. 

Consider two surfaces, one is the horizontal plane z = 0, and the other 
is defined by the equation z — (see Figure 2). Consider also two vertical 
segments PiQi and P2Q2 connecting the two surfaces. We propose the problem 
of finding the surface of minimal area, bounded by these two segments and the 
two surfaces. A surface with such characteristics is known in mathematics as a 
minimal surface.^" From the point of view of Physics, the surface corresponds 
to a soap film. This follows from the fact that the soap film must minimize the 
elastic energy, what is equivalent to minimize the area of the surface. 

We will approximate this soap film by a vertical surface E, that is, a surface 
that is projected orthogonally on a curve 7 connecting points Qi and Q2 in the 
plane Oxy. As we shall see, this approximation fits very well the experimentally 
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observed soap film.^^ 

The area of the surface E for a given curve 7 is: 




where the element of area dA is the product of the length element ds = ^1 + y'^dx 
of the curve -f by (see Figure 2). 

We observe that Eqs. ([2]) and ^ are identical except for a constant factor. 
This means that the curve that minimizes the area ^(7), is the same curve that 
minimizes the time of descent ^(7). Therefore, the solution curve is the same as 
for the Brachistochrone problem, that is, a cycloid. 

Figure 2 shows a representation of the surface S and the cycloid 7. The 
parametric equations of the cycloid in the coordinates of the figure are: 

X = xq + R{a — sin a) ; y — — cos a) : < a < 27r, (4) 

where Xq and R can be determined by the coordinates of the points Qi and Q2- 
We observe from Eq. ^ that ^(7) is also the length of the curve 7 with the 
metric given by ds* — -^ds, that is, the Euclidean metric ds ~ ^ dx^ + dy"^ 
divided by y^. The above discussion shows that the cycloid is the curve along 
which the distance between two points Qi, Q2 for the metric ds* is minimized. 
For this reason, ds* is known as the cycloidal metric. 

The perspex model (a) in the photograph of Figure 3 corresponds to Figure 2. 
It consists of two perspex sheets with the profile oi z = , which are the sides 
of the model; they are glued to a perspex square base of 21cm of side, and, 
along the curved profiles, to a flexible transparent plastic sheet. We have used 
perspex sheets of 2mm of thickness and the profiles have been made by laser 
cutting in a plastics workshop. The two transversal pins can be toothpicks or 
plastic wires. 

When our model is dipped into a soap solution, a soap film is formed between 
the two pins. Incidentally, if soap films between the pins and the lateral sides 
are also established, they should be removed by piercing with a stick. We have 
drawn on the base of our model, the cycloid that rolls over the x-axis and passing 
through the base points of the pins. We can see in the photograph of Figure 4, 
that the horizontal section at 2: = of the soap film fits very well to the drawn 
arc of cycloid. ""^^ 

In our experiment, to minimize the gravity effects, we have assumed that the 
perspex base of our model has been placed horizontally. It can be instructive to 
see how the soap film is deformed when the base of the model is placed vertically. 
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Figure 3: Perspex models. 



3 The solution to a generalized Steiner problem 
using soap films. 

The classical Steiner problem consists of finding a point Q in the plane such 
that minimizes the sum of its distance to three given points Qi, Q2, Qs, that is, 

d{Q)^d{Q,Qi) + diQ,Q2)+diQ,Q3) (5) 

is minimized, where d(Q, Qi) is the Euclidean distance between Q and Qi, {i = 
1, 2, 3).'' We assume that each angle of the triangle Q1Q2Q3 is less than 120°. 

This problem, that has many industrial and commercial applications, has 
an experimental solution using soap films. ^'^ To see this, consider the soap film 
limited by two parallel plates placed horizontally and three vertical pins. The 
soap film is composed of three vertical planes, which meet in a vertical edge 
over a point Q at angles of 120° to each other, according to the second Plateau 
law.® Because the area minimization property of this soap film is equivalent to 
length minimization property of the base curve, this soap film gives the solution 
to the Steiner problem. 

Consider now the Steiner problem for the cycloidal metric, so that we have 
to minimize the sum of the length of three cycloidal arcs that meet at certain 
point Q. The solution to this problem can be obtained experimentally by using 
soap films. For that, we only have to add to the perspex model described in 
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Figure 4: Photographs of the soap film formed between surfaces z = 0, z = 

and the vertical pins over points Qi — (2.3,2.3) and Q2 = (13.6,3.3) in cm. Its 
intersection with z — fits very well the arc of cycloid passing through these 
two points and with base line the x-axis. This cycloid is generated by a circle 
of radius R = 2.54 cm rolling on the x-axis. 
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Figure 5: Photograph of soap film solution to the Steiner problem for the cy- 
cloidal metric. Note that the three film components meet at angles of 120° to 
each other. 



Section [5] a third vertical pin. Then, the base curve of the soap film that sets 
for this model, is the solution to this Steiner problem. 

Figure 5 shows this experimental solution. Observe that the film components 
meet at angles of 120° as in the classical Steiner problem. 

4 Solution to the problem of quickest approach 
with soap films. 

The solution of Jakob Bernoulli to the Brachistochrone problem was published, 
with the other received solutions, in 1697 in Acta Eruditorum under the title 
"Solution of problems of my brother together with the proposition of other in 
turn" . 

In this article he proposed two new variational problems. The first is known 
as the problem of quickest approach, and the second is related with the isoperi- 
metric problem. 

Jakob enunciated the first problem in these terms (see Figure 6): 

On which of the infinitely many cycloids passing through O, with the same 

base OA, can a heavy body fall from O to a given vertical line AB in the shortest 

time? 

His brother Johann solved this problem quickly, showing that the cycloid 
of quickest approach is the cycloid that meets orthogonally with the vertical 
straight line AB (see Figure 6). He also solved the problem for the general case 
in which AB is not a straight line. The solution is, also in this case, the cycloid 
orthogonal to AB.^^ 

In Figure 7 we can see the solution to this problem using soap films. We 
have used the model described in Section [2] consisting of two perspex plates 
corresponding to the surfaces z — and z ~ and we have consider again, 
as a good approximation, that soap films between vertical pins or planes are 
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Figure 6: The family of cycloids passing through O and with the same base OA. 
The cycloid along which a heavy body falls from O to AB in the shortest time 
is precisely the arc OQP that meets AB orthogonally. The transversal curves 
are the synchronies, that is, the curves of the simultaneous positions of heavy 
particles which are released at along the cycloids at the same time. 



Figure 7: Photographs of soap film solution to the problem of quickest approach. 
The pin is over the point (4, 5.5) in cm. The soap film meets the vertical surface 
orthogonally (1st Plateau law), and its horizontal section at z = is an arc of 
cycloid with R = A cm. 
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Figure 8: Vertical surface S of minimal area bounded by z = 0, z = ky, and 
two vertical pins. Its horizontal section at z = is a catenary 7. 



vertical surfaces. Translated to this model, the problem is to find the surface 
of minimal area that connects a given vertical pin with a given vertical surface. 
The soap film solution corresponds approximately to a vertical surface with base 
a cycloidal arc that meets the given vertical surface at an angle of 90°. The 
configuration for other angle is not stable, as expressed by Plateau first law.* 
The required arc of cycloid (see Figures 6 and 7) has a generating circle whose 
perimeter is twice the distance OA. Note that AP = 2R, OA = ttR. 



5 The catenary and the Poincare half-plane. 

5.1 The catenary. 

In 1691 Jakob Bernoulli proposed as a challenge to find the curve that assumes a 
freely hanging chain when held down by its ends. His brother Johann, Christiaan 
Huygens, and Gottfried Leibniz found the curve solution, since then known as 
the catenary. 

As a variational problem, the catenary has to minimize the gravitational 
potential of the chain, that is: 



^(7) = / gydm = gp yds ^ gp y^/l + y''^dx (6) 

where gydm is the potential energy in an uniform gravitational field g of an in- 
finitesimal element ds with mass dm and density p, which is situated at altitude 
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Figure 9: Photographs of the soap film formed between z — 0, z — ky, and the 
vertical pins over the points (—5.6,-15.2) and (5.6,-15.2) in cm. Its intersec- 
tion with 2 = fits very well the catenary y = acosh(|) with a = 2.13 cm. This 
is the stable solution (minimum energy) . We have represented also the unstable 
solution (maximum energy) corresponding to a = 14.1 cm. 



y with respect to the point of zero potential energy level of reference. 

The explicit equation for a chain of length L, two suspension points at the 
same altitude and with separation iJ, can be written, for an adequate choice of 
the axis, as: 

(7) 



y 



X 

a cosh — 
a 



where a depends on H and L}^ 

On the other hand, the area of the vertical surface E, raised over a curve 7 
(see Figure 8) is given by 



^(7) 



ds = k y\/ 1 + y'^dx 



(8) 



Since Eqs. dU and dH) only differ in a constant factor, the curve 7 that 
minimizes 1^(7) is the same that minimizes ^(7), namely the catenary. 

If in our perspex model of Section 2 we replace the top curved surface z = 
by the plane z = ky with slope k = tan 15° (see Figure 3 (b)), the soap film 
formed between two vertical pins is approximately a vertical surface whose base 
curve is a catenary (see Figures 8 and 9).^^ There are two solutions correspond- 
ing to a minimum and a maximum of elastic energy, so the first solution is stable 
while the second is unstable.^'' 

In our case, for the vertical pins over the points (—5.6, —15.2) and (5.6, —15.2) 
in cm, the two solutions calculated from Eq. [7] correspond to a = 2.13 cm and 
a = 14.1cm. The photographs of Figure 9 show that the experimental soap film 
fits very well the catenary drawn on the bottom plate according to the above 
equation for the stable case. 

David Lovett also studied the catenary with the help of soap films by show- 
ing that the soap film formed between two pins across two non-parallel plates 
forming a wedge is part of a catenoidal surface.^' 
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Figure 10: Vertical surface S of minimal area bounded by z = 0, 2 = i, and 
two vertical pins. Its horizontal section 7 is an arc of circle with center in the 
X-axis. 



5.2 The Poincare half-plane. 

After demonstrating that the cycloid is the curve that minimizes , Johann 

Bernoulli showed that the curve 7 that minimizes f — is an arc of circle. 
This means that, considering the Poincare half-plane, that is, the half-plane 
— {{x,y),y > 0} with the metric ^ = \^dx +dy ^^ curve of minimal 
longitude joining two points Qi and Q2 of i7+ that are not in a vertical line, 
is the only semicircle passing through Qi and Q2 and with center a point C on 
the axis y = (see Figure 10).-^^ 

Again, we can get the above solution using soap films. For that, we have to 
replace the top surface of our perspex model with a surface of equation z — 
(see Figure 3(c)), so that the soap film that is formed between two vertical pins 
is approximately a vertical surface with base curve a semicircle (see Figures 10 
and 11). This follows from the fact that the area of a vertical surface over 7 is 
given by: 

7-7 2/ y 

and therefore minimizing the area of S is the same as minimizing the length of 

7 in the metric — . 

' y 

We can see in Figure 11 that the experimental soap film fits very well the 
semicircle drawn in the bottom plate. 

The soap film in the Poincare half-plane model is specially interesting for 
visualizing geometric properties of the cosmological model of the universe associ- 
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Figure 11: Photographs of the soap film formed between z = 0, z = -i, and two 
vertical pins over the points (—5.6, 3.3) and (5.6, 3.3) in cm. Its intersection with 
z = fits very well an arc of circle centered in the x-axis with radius i? = 6.5 
cm. 



ated to a matter density less than the critical, because in this case the space has 
negative constant curvature so that the Poincare half-plane is a bidimensional 
model for this space. 

Here, as in the case of the cycloid, we can consider the soap film solution 
in terms of the Fermat principle for geometric optics, and consider the arcs of 
catenary and the arcs of circle as the path of the light rays in a medium with 
refractive index proportional to y and ^ respectively. We can also study soap 
films formed by three or more vertical pins, as experimental solutions of the 
associated Steiner problems for the metrics yds and y. The last case is of 
special interest because it gives an experimental solution to the Steiner problem 
in the Poincare half-plane space. 

Finally, it is interesting to observe that the cycloids and the arcs of circle 
can also be considered catenaries for gravitational potentials proportional to 

and ^, respectively. 

6 Summary 

We have seen how to design experiments with soap films formed between two 
perspex surfaces to study the solution to several variational problems that ap- 
pear in mechanics, optics, elasticity, and geometry. 

We have studied the brachistochrone problem, the Steiner problem, the prob- 
lem of quickest approach, the catenary problem and the Poincare half-plane 
problem. All are very important from a pedagogical point of view. 

In our perspex models we have taken one of the surfaces as a plane (z = 0) 
and for the other have considered three possibilities: z = z = ky, and z = i, 
such that the soap films formed give the solution to the brachistochrone problem, 
the catenary problem, and the Poincare half-plane problem, respectively. To 
prove this we have seen that the problem of minimizing certain mechanical 
magnitudes is equivalent to minimizing the area of a surface connecting two 
pins between two appropriate surfaces. 
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The above problems can also be considered as the optical problems of finding 
the path that follows a light ray, moving in a medium whose refractive index is 
proportional to y, or 1/y, respectively. In the context of geometry these 

problems consist in to calculate the geodesies for the metrics ds/ yfy, yds, or 
ds/y, respectively, where ds = \/ dx^ + dy^ is the Euclidean metric. Finally, 
the curves solution to the above problems (that is, arcs of cycloid, catenary 
and circle), can all be considered catenaries by taking gravitational potentials 
proportional to j/, and i, respectively. 

These considerations illustrate in a simple way, the beautiful connection 
between dynamics, statics, optics, elasticity, and geometry. 

Moreover, we have pointed out the historical origin of these problems, with 
the challenges of the brothers Jakob and Johann Bernoulli. 
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Appendix A: Suggested problem 

Consider the following two surfaces: one is the horizontal plane z = 0. and 
the other is defined by z = y^, y > 0. Consider also two vertical segments 
connecting these surfaces. 

(a) Prove that the vertical surface of minimal area bounded by these two seg- 
ments and the two surfaces intersects the plane ^ = in an arc of parabola. 
For which positions of the segments there is not any solution? (Hint: The 
Euler-Lagrange equation gives two solutions but only one corresponds to 
a minimum.) 

(b) Find the refraction index for which the light rays in iJ+ = {(x, y); y > 0} 
describe parabolas. 

(c) Find the gravitational potential for which the catenary is an arc of parabola. 

(d) What is the metric in ff + for which the parabolas are the geodesies? 



°^ Electronic mail: c_criado@uma.es 

Electronic mail: alamo@uma.es 

^ Johann Bernoulli challenged the mathematical world with the 
brachistochrone problem, published in Acta Eruditorum in June, 1696, see D. 
J. Struick (ed), A source book in Mathematics, (Princeton University Press, 
Princeton, New Jersey, 1986), pp. 391-392. 
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^ To obtain this sohition, he used the Fermat Principle and considered that the 
curve of shortest descent time must be the same that the curve described by a 
light ray whose velocity is proportional to the square root of the altitude, that 
is, the velocity that a heavy body acquires in falling. This sohition was 
published in Acta Eruditorum in May, 1697, pp. 206-211. In the same issue of 
that journal, appeared also the solutions to the brachistochrone problem given 
by his brother Jakob Bernoulli, Leibniz, L'Hopital, Tschirnhaus, and Newton. 
All except L'Hopital find the cycloid as the solution. The method used by 
Jakob Bernoulli contained the roots of the calculus of variations (see reference 
in Note 1). 

^ J. A. F. Plateau, Statique experinientale et theorique des liquides soumis aux 
seules forces moleculaires, 2 Vols. (Gauthier-Villars, Paris, 1873). 

* R. Courant and H. Robbins, What is Mathematics? (Oxford University 
Press, London, 1941) 

^ C. V. Boys, Soap-Bubbles, their colours and the forces that mould them 
(Dover, New York, 1959). 

^ S. Hildebrant and A. Tromba, Mathematics and Optimal forms (Scientific 
American Books, New York, 1985). 

^ C. Isenberg, The Science of Soap Films and Soap Bubbles (Dover, New York, 
1992). 

^ D. Lovett, Demonstrating Science with Soap Films (lOP-publishing, Bristol, 
1994). 

^ A sufficient condition for a curve 7 to make the integral /(y, y' , x)dx 
extremal is that it satisfies the second order differential equation 
1^ — = 0, known as the Euler-Lagrange equation (see H. Goldstein, 

Classical Mechanics (Addison- Wesley, Reading, Massachusetts, 1980), p.45. 
For the case in which / does not depends on x, the Euler-Lagrange equation 

reduces to y' §^ — f = constant. In our case, /(y, y') = and then we 

get the differential equation y' = \J^-^t where fc is a constant. It is easy to 
see that its sohition is the cycloid. 

R. Osserman, A Survey on Minimal Surfaces (Dover Publications, New 
York, 1986). 

According to the first Plateau law (see, for example, Ref. 8. pp. 8-9), the 
film have to meet both surfaces z = Q and 2; = ^ at an angle of 90°. 
Therefore, the approximation of taking a vertical surface S is better as y is 
larger, because then the slope of the tangent to 2: = 1/ ,Jy is smaller. 

^"^ A rough estimation of the relative error when we take a vertical clement of 
surface instead of the element of surface for the soap film can be 
calculated by 

AA' - AA (a - tan a) As 

= IT-, = (10) 

A^ tana As ^ ' 
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where AA' is the element of surface for a cylindrical surface that meets 
orthogonaly the two plates, As is the line element along the curve 7, and 
tana = ^ = z' . 

Neglecting the terms of order greater than 4 in the expansion of a — tan a as 
function of the derivative z' , we get — — which for z = 1/ ^ leads to 
Cr — —y~^/12. Therefore, for y > 2 we have an error less than 1%. 
Similar calculations can be made for z = ky and z = l/y for the cases of the 
catenary and the Poincare half-plane, respectively. 

C. Isenberg, "Problem solving with soap films," Phys. Educ. 10, 452-456 
(1975). 

^'^ See reference in Note 1, p. 399 

Johann Bernoulli could easily solve this problem because of his early study 
on the synchronies, which are the orthogonal trajectories to the cycloids 
passing through O (see Figure 6). From a mechanical point of view, these 
curves are formed by simultaneous positions of heavy particles which are 
released at O along the cycloids at the same time. On the other hand, from an 
optical point of view, the synchronies correspond to the simultaneous positions 
of the light pulses emitted from O at the same instant and that propagate 
along a medium with refractive index proportional to 1/ y^. Therefore the 
synchronies are the wavefronts. From the fact that wavefronts and light rays 
are perpendicular, Johann Bernoulli concluded that the synchronies are the 
orthogonal trajectories to the cycloids. Thus the point of quickest approach P 
must be the point where one of the synchronies is tangent to the line AB (see 
Figure 6). For more details on this story see S. B. Engelsman, Families of 
Curves and the Origins of Partial Differentiation (North- Holland, Amsterdam, 
1984), pp. 31-37. 



The parameter a is determined by ^ + {^)^ — a cosh The suspension 
points are symmetric with respect to the y-axis and their distance to the x-axis 
is & = ^ + (-j)^. For the proof of these facts see J. L. Troutman, Variational 
Calculus and Optimal Control (Springer Verlag, New York, 1996), pp. 78-80. 

When the ratio between pin separation and distance from the pins to the 
X-axis, is greater than 1.335 there is not any parameter a verifying 
6 = a cosh 1^ (see Note 16), while for less than 1.335 there are two values of 
parameter a. 

^® The catenoid is the surface of revolution generated by the catenary, and, as 
Euler showed in 1744, is the only minimal surface of revolution (see Ref. 6). 
The reason is that the area of a surface of revolution generated by a curve 7 of 
equation y — y{x) is given by ^(7) = 27r yds = 2tt y\/l + y'^dx, and 
again coincides, except for a constant factor, with Eqs. ([6]) and dS]). The 
Euler-Lagrange equation (see Note 9), reduces in this case to = 

and its integration gives the catenary of equation y — ci cosh . 
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^ The Euler-Lagrangc equation (sec Note 9), reduces in this case to 
— . ^ = Ci , and its integration gives the semi-circle of equation 

y = ci sin t,x = C2 + ci cos t, for < t < tt. 

The Poincare half-plane is one of the simplest models of non-euclidean spaces. 

The geodesies of this space are the semicircles with center a point C on the 
axis y = together with the vertical straight lines. If we think of these 
geodesies as the "straight lines" of this space, the associated geometry does 
not verify the 5*'* Euclidean postulate (the parallel postulate), because given a 
"straight line" and a point not on that line, there are infinitelly many lines 
through the point that do not cut (i.e. parallels to) the given straight line. 
Therefore, only a small step, considering semicircles as "straight lines" 
separated Johann Bernoulli to discover non Euclidean geometry 130 years 
before Bolyay and Lobatchesky did. 

W. Rindler, Relativity: Special, General and Cosmological, (Oxford 
University Press, New York, 2001), p. 364. 
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